Quantum phase transitions in the extended periodic Anderson model 



Akihisa Koga, 1,2 Norio Kawakami, 1 Robert Peters 2 and Thomas Pruschke 2 
1 Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
Institut fur Theoretische Physik Universitdt Gottingen, Gottingen D-37077, Germany 



(Dated: February 1, 2008) 

We investigate quantum phase transitions in the extended periodic Anderson model, which in- 
cludes electron correlations within and between itinerant and localized bands. We calculate zero and 
finite temperature properties of the system using the combination of dynamical mean-field theory 
and the numerical renormalization group. At half filling, a phase transition between a Mott insulat- 
ing state and a Kondo insulating state occurs in the strong coupling regime. We furthermore find 
that a metallic state is stabilized in the weak coupling regime. This state should be adiabatically 
connected to the orbital selective Mott state with one orbital localized and the other itinerant. The 
effect of hole doping is also addressed. 



I. INTRODUCTION 

Strongly correlated electron systems with orbital de- 
generacy have attracted considerable interest. One typi- 
cal example is the transition metal oxide L1V2O4. In this 
compound, unexpected heavy fermion-like behavior was 
observed^ which is quite likely induced by an interplay 
between the geometrical frustration of the spinel struc- 
ture and the orbital degeneracy of the d electrons 
One particular scenario suggests that due to the V posi- 
tions in the local trigonal crystal environment, the ti g 
states at the Fermi energy are split into a\ g and 
subbands, which then play different roles in stabilizing 
the heavy electron state.— ^ For example, due to different 
bandwidths of the a ig and e g bands, the local Coulomb 
correlations can induce a Mott insulating state in the 
former while the latter stay metallic^ A hybridiza- 
tion between the two subshells on neighboring V atoms 
in the unit cell then can explain the observed heavy 
fermion behavior. Other examples are the compounds 
Ca 2 - x Sr x Ru04— and La n +iNi n 03 n +iriii^ where a sim- 
ilar orbital-selective Mott transition-^ is suggested to be 
induced by the chemical substitutions and the change of 
temperature, respectively. It was furthermore reported 
that a heavy fermion state is realized in the former com- 
pound around x « 0.2^2- while a bad metal with localized 
spins at low temperatures is characteristic of the latter ^ 
These observations stimulate further theoretical investi- 
gations on multiorbital system o 14 ' 15 ' 16 ' 17 ' 18 ' 19 

A common feature of those interesting examples is that 
localized and itinerant electrons and their correlations in 
the multiorbital system seem to play an important role in 
the formation of heavy fermion states in transition metal 
oxides. Motivated by this, we here investigate a multi- 
orbital system with localized and itinerant bands. It is 
known that the hybridization together with local elec- 
tron correlations leads to the Kondo effect and in turn 
to the large density of states at the Fermi level, which 
explains the heavy fermion states observed in rare-earth 
compounds ! 20 ' 21 ! 22 On the other hand, electron correla- 
tions in an itinerant band also yield a heavy mass and 



eventually induce a transition to the Mott insulating 
state d 3 ' 28 Thus the interesting question arises whether 
the heavy fermion states originating from the Kondo ef- 
fect and close to the Mott transition can be distinguished 
from each other or not. This question has recently been 
discussed by several authors ; 14 ! 15 ' 24 ' 25 but it was not clear 
how the heavy fermion state is realized at low tempera- 
tures. In particular, the role of the Hund coupling needs 
to be clarified, which may be important in /-electron sys- 
tems as well as transition metal oxides. Therefore, it is 
necessary to discuss low temperature properties in cor- 
related multiorbital electron systems with localized and 
itinerant bands systematically. 

To this end, we investigate an extended version of the 
periodic Anderson model, where not only onsite Coulomb 
interactions in the localized and itinerant bands but also 
the interband Coulomb interaction and Hund coupling 
are taken into account. The properties we are inter- 
ested in - heavy fermion behavior and Mott transition 
- are well captured by the dynamical mean field the- 
ory (DMFT)^! 2 !^^ and furthermore do not depend 
critically on details of the band structure. Thus we will 
discuss the low-temperature properties for a particle-hole 
symmetric system within DMFT. We will clarify that the 
Hund coupling between different orbitals plays an impor- 
tant role, leading to a quantum phase transition between 
a Kondo insulating phase, a Mott insulating phase and 
the metallic phase. The Mott and Kondo insulators also 
behave differently when doping the system. 

The paper is organized as follows. In Sec. UH we intro- 
duce the extended periodic Anderson model and briefly 
summarize DMFT and numerical techniques. The com- 
petition between some phases is discussed in Sec. Mil 
We also discuss the effect of hole doping in Sec. US A 
brief summary is given in Sec. [V] 

II. MODEL AND METHOD 

We consider a correlated electron system with orbital 
degeneracy described by the Hamiltonian 
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where c| Q0 .(ci Q , cr ) creates (annihilates) an electron with 
spin (r(=t, J,) and orbital index a(= 1, 2) at the ith site, 
and n iaa — c\ aa c iaa . For orbital a, represents the 
transfer integral, V the hybridization between orbitals, 
U a and U' the intra-orbital and inter-orbital Coulomb 
interactions, J the Hund coupling, and \x the chemical 
potential. The structure of the model Hamiltonian is 
schematically shown in Fig. [T] We note that the model is 
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FIG. 1: Sketch of the structure of the model Hamiltonian 

reduced to conventional systems in several limiting cases. 
For example, when V = 0, the system becomes the two 
orbital Hubbard model, where we expect a Mott tran- 

(2) 

sition to occur. On the other hand, when iV. = 0, we 
recover at U\ = U' = J = the conventional periodic 
Anderson model for heavy-fermion systems^ 21 - while the 
choice Ui = U2 = V = 0, J > leads to the double 
exchange model discussed in connection with magnetism 
in transition metal oxides i 30 i 31 ' 32 i 33 

In this paper, we want to focus on the low-temperature 
properties of a system with one band localized and 

one kept itinerant, i.e. we choose the hopping integral 

(2) 

for the a = 2 band as th = 0, but allow all other 
parameters to be finite. To treat the extended peri- 
odic Anderson model in this parameter regime, we use 
the DMFT. This method has been developed in sev- 
eral group o 26 i 27 ' 28 i 29 and has successfully been applied 
to correlated electron systems such as the single-band 
Hubbard model^^^i 3 !^^^ the two-band Hubbard 
mode&l^L^Md^je^s Qr the periodic Anderson 

model^i 2 ^^^^ 2 .^ 3 . In the DMFT, a lattice model is 
mapped to an effective impurity model, where local elec- 
tron correlations are taken into account precisely. The 
requirement that the site-diagonal lattice Green function 
is equal to the one of the effective quantum impurity 



then leads to a self-consistency condition for the param- 
eters entering the impurity problem. This treatment is 
formally exact for infinite spatial dimensions. If one is 
allowed to ignore nonlocal correlations, e.g. sufficiently 
far away from phase transitions with long-range order, 
the method can be used as a quite reliable and accurate 
approximation to three dimensional systems. 

Within the DMFT, the lattice typically enters only 
via the density of states (DOS) of the system with 
vanishing two-particle interactions. Since we are inter- 
ested in generic features of the extended periodic An- 
derson model, we are free to choose a convenient lat- 
tice structure. Here, we use a Bethe lattice with infi- 
nite coordination^ for which the DMFT self-consistency 
equation is simplified to 24 ' 25 



G toc {z) 
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where Go is the non-interacting Green function for the 
effective impurity model and Gi oc is the local Green func- 
tion. Both are matrices with respect to the orbital index, 
but due to z-j = only the (1,1) element is needed. 
As is well-known^ the relation (0| is equivalent to a 
semi-elliptic DOS for the itinerant band (a = 1) with 
half-bandwidth D. 

Let us note that the localized band (a — 2) does not 
appear explicitly in the self-consistency condition eq. ([J). 
For the effective impurity model entering the DMFT, this 
means that only one of the impurity orbitals connects to 
the rest of the system through an effective hybridization 



function defined by 
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In the following, 



Gi oc (z) 

this quantity is represented by an auxiliary set of dynam- 
ical degrees of freedom, which we will call an effective 
bath. 

The Hamiltonian representing the effective impurity 
model can now be written as 



k ka 

+ Hi 

ocal 1 (5) 



where a ka are the auxiliary annihilation operators with 
quantum number k and spin a defining the effective bath, 



3 



E a = —fjt is the energy level for the impurity site, n aa — 
c qct c q(t and H[ oca i the local interaction term as defined in 
(J3j) . The set of parameters for the effective bath {ek,Jk} 
must be determined self-consistently through the DMFT 
condition eq. (|4]). 

The Hamiltonian ([5]) represents a quantum impurity 
model, which is a challenging theoretical problem on its 
own account. Thus, to discuss low-temperature prop- 
erties of the extended periodic Anderson model in the 
DMFT, we need a tool to accurately solve such quan- 
tum impurity models. While applicable in the weak- 
coupling regime, perturbation theory in general fails in 
the strong coupling regime, where for example the antici- 
pated transition between Mott and Kondo physics should 
appear. Exact diagonalizations and quantum Monte 
Carlo (QMC) simulations are known to be numerically 
exact methods, but both cannot properly resolve expo- 
nentially small energy scales, which again will play an im- 
portant role in the interesting regime, where Kondo and 
Mott physics compete. In addition, due to the required 
computational resources, QMC simulations typically can 
access only a restricted window of the parameter space 
and can also be subject to a severe sign problem when 
applied to complex multi-orbital models. 

Therefore, we here use Wilson's numerical renormal- 
ization group^ (NRG) to solve the effective impurity 
model. In the NRG, one discretizes the effective bath on a 
logarithmic mesh by introducing a discretization param- 
eter A > 1. The resulting discrete system can be mapped 
to a semi-infinite chain with exponentially decreasing 
couplings^ which allow us to access and discuss prop- 
erties involving exponentially small energy scales quanti- 
tatively. To ensure that sum rules for dynamical quanti- 
ties are fulfilled, we use the complete-basis-set algorithm 
proposed recentlyj 55 i 56 ' 57 We observe that obeying the 
sum rules is mandatory to properly describe the low en- 
ergy properties of the system away from half filling. In 
the NRG calculations, we use a discretization parameter 
A = 2 and keep 1200 states at each step. 



III. RESULTS AT HALF FILLING 

In the following, we will present results calculated for 
U = Ux = U 2 with J/U = 0.1 and U' = U - 2 J fixed. 
The effect of Hund's coupling will be addressed later. 

Let us begin with the low-temperature properties of 
the extended periodic Anderson model at half filling, 
which is realized by fixing the chemical potential fi = 
/Jo = U/2 + U' — J/2. For the conventional periodic An- 
derson model Ux = U' = J = 0, it is known that the 
band insulator at (V/D ^S> 1) is adiabatically connected 
to the Kondo insulator with U2/D 3> 1. Furthermore, it 
was claimed that the introduction of a Coulomb interac- 
tion for the conduction band U\ does not induce a tran- 
sition to the Mott insulating state i 24 ' 25 Namely, in this 
case always antiferromagnetic spin correlations between 
the orbitals develop. However, in the presence of Hund's 



coupling J, which induces ferromagnetic spin correlations 
between the orbitals, one should expect a competition be- 
tween the singlet formation due to the hybridization and 
a high spin formation due to the Hund coupling. 

To discuss how spin correlations develop between the 
degenerate orbitals, we calculate the squared spin mo- 
ment (S 2 ), where S z = ( n a] — n ai) /2, as shown in 
Fig. [2 When V/D is large, the local spin moment is 
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FIG. 2: (color online) The local squared moment (S z ) as a 
function of the hybridization V/D. Solid (open) circles repre- 
sent the results at temperature T/D — 2.1xl0~ 2 (6.5 xlO -4 ). 

strongly suppressed, (Si) — ► 0. In this regime, the hy- 
bridization together with the intraorbital Coulomb inter- 
actions leads to the formation of an interorbital singlet, 
where the insulator (Kondo insulator) will be realized in 
the half-filled particle-hole symmetric system. On the 
other hand, different behavior occurs for small V/D, as 
shown in Fig. [2j Here we find that (S" 2 ) — > 2/3 when 
U/D > 2, while (S 2 ) -> 1/2 in the other limit. 

In the former case, strong correlations localize the elec- 
trons in band a = 1. Note that due to a finite Hund's 
coupling, the corresponding Mott transition will occur 
at a smaller value of U as compared to the one-band 
Hubbard model here . 1 Furthermore, for small V, the 
Hund's coupling will win over the effective antiferromag- 
netic exchange ~ V 2 /U generated by the hybridization 
and the Coulomb interaction, i.e. a spin triplet state 
S = 1 is realized, which implies (S 2 ) =2/3. Since for 
large V/D the ground state is a local singlet with 5 = 0, 
we expect that at T — a phase transition from the Mott 
insulating state to the Kondo insulating state occurs as 
a function of V/D. In fact, the phase transition is clearly 
visible for the lower temperature, as shown in Fig. [2] 

On the other hand, in the regime U/D < 2, we ob- 
serve a smooth curve for (S 2 ) as a function of V, which 
furthermore depends only little on temperature, except 
for a small region around the turning point. This be- 
havior is shown in Fig. [5] for U/D = 1 as example. For 
small V/D Hund's exchange will again dominate over the 
antiferromagnetic interorbital coupling generated by V. 
Thus we are effectively left with a Fermi liquid ferromag- 
netically coupled to a localized spin, i.e. a ferromagnetic 
Kondo model. As is well-known, this model has a metal- 
lic ground state,™ where the local spin acts as a potential 
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scatterer. The increase of V eventually leads to a domi- 
nance of the hybridization term and a Kondo screening of 
the local spin by the itinerant electrons. For finite tem- 
peratures one will observe a crossover, because only for 
Tk > T a full singlet with (5 2 ) = between local and 
band spins is realized. With decreasing temperature, the 
crossover becomes sharper and, because one now probes 
lower Kondo temperatures, shifts to the left, as seen in 
Fig. [2j Eventually it turns into a transition at T = 0, 
where the Kondo state always leads to a singlet for arbi- 
trary small antiferromagnetic coupling. 

This qualitative discussion can be further substanti- 
ated by looking at the density of states for itinerant and 
localized orbitals in Fig. [3] We first focus on the case 
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FIG. 3: (color online) Density of states for the extended 
periodic Anderson model with fixed U/D = 1.0 (a) and 
U/D = 3.0 (b) at T/D = 6.5 x 10" 4 . The numbers repre- 
sent the values V/D. 

U/D = 1. As anticipated before, for small V/D a metal- 
lic state is realized in the orbital a = 1. The orbital 
a = 2, on the other hand remains localized. To confirm 
the stability of the metallic state and support the qual- 
itative discussion, we also show the density of states at 
a much lower temperature T/D ps 10~ 8 as the dashed 
lines in Fig. [5] As expected, we find that the metal- 
lic state is certainly the ground state for the a = 1 or- 
bital, while a strongly suppressed DOS at Fermi level is 
seen in the other orbital. The NRG results are consis- 
tent with those obtained from a calculation using exact 
diagonalization ; 28 ! 34 see full lines in Fig. [4] Further in- 
crease of the hybridization beyond the phase boundary 
(V/D) c ps 0.14 eventually leads to the appearance of a 
Kondo resonance in P2{^), as shown in Fig. [3] (a), and a 
corresponding reduction of the DOS at the Fermi energy 
in pi(pS). For even larger V, one then finds the structures 
known from Kondo insulatorS ) 49 i 52 ' 58 i.e. a Kondo reso- 
nance split by a hybridization gap in p2(w) and a gap of 




FIG. 4: (color online) Solid and dashed lines represent the 
density of states obtained by exact diagonalization with 12 
sites and the NRG method at T/D » 10" 8 for the metallic 
state with U/D = 1.0 and V/D = 0.1. 



the same width in pi(ui). Note that one would actually 
expect a similar gap to appear in p^i^j) at V/D = 0.16. 
However, this gap appears only for T -C Tk and is thus 
not yet visible here. 

When U/D = 3, on the other hand, we always find an 
insulating state, which can be characterized as a Mott 
(Kondo) insulator with a charge gap in both orbitals for 
small (large) V/D. When V approaches the critical value 
(V/D) c w 0.42 from below, the charge gap around the 
Fermi level decreases, but does not seem to vanish even 
at the critical point, which implies that the phase tran- 
sition is of first order, in accordance with the behavior 
of (5 2 ). The initial reduction of the charge gap with in- 
creasing V can be interpreted in the following way. For 
small V, the Mott state is stabilized at a lower critical 
U c due to ferromagnetic Hund's coupling ] 43 ! 48 With in- 
creasing V, the additionally generated antiferromagnetic 
exchange effectively reduces this effect, i.e. the orbital 
a = 1 is driven closer to its critical value and the charge 
gap decreases. As soon as the antiferromagnetic coupling 
dominates, the gap scale is set by the Kondo scale. The 
properties very close and at the critical point are thus 
determined by this subtle competition of energy scales. 
Whether the two regimes are separated by a quantum 
critical point or the transition is driven by a simple level 
crossing remains to be investigated. 

By performing similar calculations for various parame- 
ter values, we obtain the phase diagram shown in Fig. [S] 
When V/D is large, the Kondo insulator prevails, where 
the spins for the localized band (a = 2) are screened by 
forming a singlet with the electrons in orbital a = 1 . The 
increase of the Coulomb interaction U with fixed J/U de- 
creases the hybridization gap characteristic of the Kondo 
insulating state. Finally, a phase transition occurs to the 
Mott insulating state with a local moment 5 = 1. Note 
that in this regime the Kondo interaction is dominated 
by ferromagnetic Hund's exchange. The competition be- 
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FIG. 5: (color online) Finite temperature phase diagram for 
the extended periodic Anderson model with J/U = 0.1 and 
T/D = 6.5 x 1(T 4 . Solid circles, open circles and solid squares 
denote the Kondo insulator, the Mott insulator and the metal- 
lic state, respectively. The dashed line is obtained from the 
simplified model (see text). 
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FIG. 6: (color online) (a) Local moment as a function of the 
Coulomb interaction U for V/D = 0.1 and 0.2. Solid (open) 
circles are the results for T/D = 2.1 x 10~ 2 (6.5 x 10" 4 ). (b) 
Density of states for the orbital a — 1 and different values 
U/D = 0.4, 0.8, 1.2, 1.6, 2.0 and 2.4 for T/D = 6.5 x 10~ 4 . 



tween these phases can be qualitatively described by a 
simplified local Hamiltonian for two orbitals, since the 
bandwidth has little effect on the phase transition. By 
examining the lowest state for this simplified model, we 
find a level crossing at V/D = \[23 '/D between the sin- 
glet and the triplet states, which is shown as the dashed 
line in Fig. [5l Obviously, this line is in good agreement 
with the phase boundary obtained by DMFT. This result 
also suggests that the actual phase transition is of first 
order and not accompanied by a quantum critical point. 

When the hybridization is small (V/D < \/2JD) and 
U/D < 2, a metallic state appears with orbital a = 2 
almost localized and a — 1 itinerant. As the Coulomb 
interaction is increased for the metallic state, the quasi- 
particle peak becomes sharper for the itinerant orbital, 
as shown in Fig. [5] Finally the metal-insulator tran- 
sition occurs to the Mott insulating state. The ob- 
tained phase boundary is consistent with the critical 
point (U/D, V/D) w (2,0) for the two-orbital Hubbard 
model with the large difference of the bandwidths. 19 

Our results are in clear contrast to those for the con- 
ventional periodic Anderson model, where at half filling 
and particle-hole symmetry, no Mott transition occurs 
even in the large U% and U2 case, but one always finds a 
Kondo insulator. Since this model is obtained in the limit 
Hund's exchange J = in our extended periodic Ander- 
son model, we can therefore conclude that the competi- 
tion between metallic and the insulating states originates 
from the Hund coupling. This observation may be rele- 
vant for the heavy fermion behavior observed in transi- 
tion metal oxides such as Ca2- x Sr x Ru04 (0.2 < x < 0.5) 
and La n+ iNi n 03„+i, which is now under investigation. 



IV. EFFECT OF HOLE DOPING 

In this section, we discuss the effect of hole doping on 
the insulating states to clarify how the heavy fermion be- 
havior emerges. To this end, we introduce the parameter 
A/i = fiQ — /i. In Fig. we show occupancy, local mo- 
ment and density of states. For (U/D, V/D) = (3.0, 0.3), 
the system is in the Mott insulating state at half filling 
A/i = 0, and the increase of A/i has little effect on the 
nature of the phase for 0.0 < Afi/D < 0.5, as shown in 
Fig. [7] (a). Further increase of the chemical potential 
eventually drives the system into a metallic state, where 
the local moment and the total occupancy n tot {— 711+712) 
both decrease. Note that most of the holes are doped in 
the a = 1 band while only a small amount of holes pop- 
ulates the other when 0.5 < A/j/Z? < 1.0. This implies 
that the character of original orbitals still survives in the 
system with hybridization. This orbital-selective behav- 
ior is also found clearly in the density of states, as shown 
in Fig. [3(b). When A/i/D > 1.3, the behavior changes: 
the number of electrons for the a — 2 band deviates from 
the half filling, as shown in Fig. [7] (a). Furthermore, we 
find that a sharp peak appears in the density of states 
characteristic of the doped Kondo insulator. This sug- 
gests that a crossover between different metallic states 
occurs around Afi/D » 1.3. 

On the other hand, when holes are doped into the 
Kondo insulating state at (U/D, V/D) = (2.0,0.4), 
monotonic behavior appears in the quantities in Fig. [7] 
(c). In this conventional heavy fermion state, 

studied in detail in the conventional periodic Anderson 
model, is realized, as shown in Fig. [7] (d). 

Before closing this section, we briefly mention the Mott 
transitions at quarter filling, which should be induced by 
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FIG. 7: (color online) Occupancy (n a ) and squared moment 
(Sf) as a function of the chemical potential (density of states) 
for V/D = 0.2 (a) [(b)] and V/D = 0.4 (c) [(d)]. 



the intraorbital and the interorbital interactions. The re- 
sults obtained for several choices of parameters are shown 
in Fig. [5] It is found that the ground state at quarter 
filling n t ot — 0.5 strongly depends on the strength of the 
Coulomb interaction U . When U is small, no remarkable 
features are visible in the curve of the total occupancy, as 
shown in Figs. [S](a) and (b). Therefore, we infer that the 
system is a metal state in the case. In the large U case, on 
the other hand, a plateau appears in the curve at quarter 
filling, as shown in Figs. 0(c), (d) and (e). In addition, 
one observes a local moment (S%) — 1/4 in this region, 
which suggests the existence of a Mott insulating state. 
Our results are similar to those for a two-orbital system 
without hybridization ) 47 ' 60 which is a Mott insulator at 
quarter filling for large enough interaction parameters. 
We can therefore conclude, that the hybridization has 
little effects on the phase diagram for the quarter filled 
system, as shown in Fig. [S](f), in contrast to that for the 
half- filled system (Fig. [5|. 
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FIG. 8: (color online) (a)-(e) The total number of electrons 
n t ot and local moment (S^) as a function of Afi/D with fixed 
parameters at T/D — 6.5 x 10~ 4 . (f) Finite temperature 
phase diagram for the quarter filling at T/D — 6.5 x 10~ 4 . 
Solid (open) circles denote the metallic (Mott insulating) 
state. 



It has been clarified that the metallic and Mott insulat- 
ing states are stabilized by the Hund exchange coupling 
between the localized and itinerant bands, which may 
be relevant to real materials, e.g. transition metal oxides 
that show heavy-fermion behavior under certain condi- 
tions. 

Although we have restricted our discussion to the para- 
magnetic state in the paper, it is an interesting and in- 
structive problem to discuss possible instabilities toward 
ordered states such as magnetic order, superconductivity, 
etc. In particular, an antiferromagnetic state induced by 
magnetic fiel d 53 ' 59 can be of considerable interest in con- 
nection with various heavy fermion compounds. These 
problems are currently under investigation. 
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